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Abstract 



The processes with three or more charged particles in the final state exhibit 

^ \ particular threshold behavior, as inferred by the famous Wannier law for (2e + 

' 5a ' ion) system. We formulate a general solution which determines the threshold 

><■ 

i-£h I behavior of the cross section for multiple fragmentation. Applications to 

several systems of particular importance with three, four and five leptons 
(electrons and positrons) in the field of charged core; and two pairs of identical 
particles with opposite charges are presented. New threshold exponents for 
these systems are predicted, while some previously suggested threshold laws 

00 ' 

' are revised. 
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I. INTRODUCTION 



The famous Wannier |1J threshold law has quite an unusual status among other threshold 
laws in quantum mechanics. Being based on an appealing mechanism, it has inspired a large 
number of studies where the law was rederived, extended, tested or rebutted. The intensity 
of these studies does not show a decrease with time as testify some representative references 
to the recent publications § § [| § @ § @ § |TJ [0 [0] (T§ |7] 
|B| f2H|j |PH| . A complete bibliography on the subject would be immense. 

In this paper we suggest a method which generalizes the Wannier mechanism when 
break-up of the quantum system on a large number (four or more) of charged fragments is 
concerned. Apparently, for the first time the particular case of the problem was treated in 
1976 in the important paper by Klar and Schlecht |22| where the threshold law was derived for 
escape of three electrons from the charged core. It was suggested that the receding electrons 
form a symmetrical configuration of equilateral triangle with the positively charged core in 
its center. The treatment was quite involving and specialized being based on hyperspherical 



coordinate system | 23| . Later Grujic [23] rederived the same result using the standard 
Cartesian coordinates where the symmetry considerations are easy to apply explicitly in full 
extent. Grujic considered also some other systems along the same lines ]25[ |26[ |27| (see 
more details in Sec.[V]). The threshold law for the three-electron escape seems to find support 
in the experimental data on the near-threshold double ionization of atoms by electron impact 
2q| . Later Feagin and Filipczyk claimed an existence of a complementary law which is 



manifested at energies somewhat above a threshold, see critical discussion in Section [V]A. 
The interest to the problem was renewed recently when two electrons and positron re- 



ceding from the core with Z = 1 charge were considered by Poelstra et al | 30|| . A brief 



note by Stevens and Feagin [[31] on complete fragmentation of H2 molecule is also to be 
mentioned. The final state in the reaction with a positron could be produced by double 
ionization of a negative ion by positron impact. However, the forthcoming experiments by 
Knudsen and co-workers |32[ concern positron impact double ionization of neutral atoms 
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[Z = 2) where the threshold law has not yet been available. Its derivation was one of 
motivations for the present study. Eventually it has developed into a general approach to 
the multi-fragmentation problem which possesses two important advantages. Firstly, our 
method describes a general situation with arbitrary number of charged fragments in simple 
terms in an arbitrary coordinate frame. Secondly, it is convenient and reliable for practical 
realizations. This allows us to clarify important conceptual aspects of the problem which 
were misunderstood or misinterpreted previously. Comparing our solution with known in 
literature results we reproduce a number of threshold exponents for different systems. At the 
same time we find that several published previously results need improvement, in particular, 
we revise the threshold law for 2e~ + e + escape. A general nature of the developed method 
is illustrated by consideration of new complicated situations with up to six charged particles 
in the final state where a number of new threshold exponents is predicted. 

In Section II we introduce particular configurations which will be called scaling con- 
figurations. They describe a multidimensional dynamic potential saddle, generalizing the 
Wannier ridge which is well known for (2e + ion) system. These configurations are related 
to rectilinear trajectories of all particles in the system and play a crucial role for the complete 
fragmentation process close to its threshold. SCs embrace the essence of previous treatments 
of particular systems, but avoid attachment to some special theoretical formalism and re- 
lated technical complications. The closest analogue of our general approach in particular 
case of three-body Coulomb systems could be found in papers by Simonovic and Grujic |33] . 



Description of small deviations from SC is given in terms of a set of harmonic oscillators 
and inverted oscillators (Section p| ). The later ones describe unstable modes which govern 
the threshold law. They are quantized following a general scheme suggested by Kazansky 
and Ostrovsky [|34] [35|| . This allows us to construct a reduced form of the wave function 
for the system of charged particles and derive the threshold law in Section IV generaliz- 
ing a procedure used previously by Kazansky and Ostrovsky [|j] [[12] for derivation of the 
conventional Wannier law. Application of a developed general scheme to some particular 
systems (Section V) is followed by concluding discussion of special features of the Wannier- 



type threshold laws, with an emphasis on a relation between the underlying statistical and 
dynamical aspects of the problem (Section VI). 

II. SCALING EXPANSION 

Our goal is to consider some atomic process which breaks an atomic particle into several 
charged fragments for low excess energy E. In this situation the motion of the fragments in 
the final state of the reaction can be described in the semiclassical approximation because 
a typical variation of the Coulomb potential [/ c ~ 1/r on the wavelengths of the fragments 
A ~ l/y/ME, 5U ~ A/r 2 , is less than a typical kinetic energy T ~ E 

5U < T 

inside the Coulomb zone r < r c = l/E where the major events take place. Therefore the 
first thing to do is to find classical trajectories which lead to the desired final state with 
total fragmentation. 

It is very important that for low energy E there exists a severe restriction onto these 
trajectories. To see this let us imagine what is happening with distances separating fragments 
when they move out of the reaction domain. If a distance separating some pair of two 
attracting fragments diminishes with time, then one should expect that this pair of fragments 
can be considered as a dipole which interacts with the rest fragments. This interaction can 
transfer the kinetic energy of the two fragments to the other fragments. Therefore one 
has to expect that eventually these two fragments will loose enough energy and form a 
bound state. If this event happens then the desired total fragmentation is not achieved. 
This discussion shows that one should look for those trajectories which exhibit a monotonic 
increase of distances separating the fragments. The point is that the lower is the available 
above-threshold energy the more restrictive this condition is. 

It is convenient to present the discussed situation considering the potential energy in the 
multidimensional configuration space where its behavior can be described as kind of "valleys 
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separated by ridges". This physical picture first suggested by Wannier for a particular 



class of reactions was discussed by Fano p6| , |37| in general case. If a system occupies some 



place on some ridge then its trajectory can either go down into some valley where a bound 
state of some fragments is created, or continue to propagate along the ridge. For the total 
fragmentation one should find a classical trajectory along a top of some ridge which leads 
from the region of small separation of fragments into the final state with infinite separation. 
It is clear that the lower is the above-threshold energy, the closer a trajectory should be to 
the top of the ridge. 

Generally speaking there might exist several such ridges which lead to the final state 
with total fragmentation. In this work we study a particular ridge, which will be called 
the scaling configuration (SC). For all the systems considered up to now we have found 
that this configuration exists. More than that, for a given system there may exist several 
different SCs. A number of examples demonstrating this property is considered in Section 
|V[ Therefore one might suspect that SC is a general feature, though this latter statement 
is not verified so far. For two electrons in the field of an ion this ridge coincides with the 
Wannier ridge. 

The basic idea is simple. As was said above one has to be sure that distances separating 
fragments monotonically increase with time. This condition is definitely satisfied if a tra- 
jectory describing N particles which have masses mj (j = 1,2, -• ■ N) obeys the following 
conditions 

Tj{t) = <f>{t)p s , j = l,2,---iV, (1) 

which are valid in the center of mass reference frame Y^j m j^j(t) = 0. We shall refer to a 
trajectory satisfying Eq.fll]) as SC. The time-independent vectors p} describe the shape of 
the SC, while the function <fi(t) gives the overall scaling factor. We will see below that this 
function monotonically increases in time thus ensuring that all distances increase as well. 
Therefore this type of motion definitely results in total fragmentation avoiding traps into 
potential valleys. It is convenient to normalize the scaling function to unity for some initial 



moment of time to 

<j>(t ) = 1. (2) 

For this normalization the vectors pj play a role of coordinates of the particles at this initial 
moment of time r}(to) — Pj i n the center of mass reference frame 

J2 m jPj = °- ( 3 ) 

3 

Notice that in SC all degrees of freedom except the only one describing the overall scaling 
factor are frozen. In this sense the SC describes a quasiequilibrium of the system. 

It is obvious that to satisfy (|I|) one should choose appropriately the initial coordinates 
Pj. Let us formulate restrictions on them. Notice first of all that in the SC the accelerations 
of the particles are 

d 2 fj(t) _ z (A , 

We presume purely Coulomb interaction, or consider Coulomb asymptote in more compli- 
cated cases, which is possible because the important distances are large (r ~ rc = 1/E) for 
low above-threshold energy. Therefore the potential energy of the system of N fragments is 

U=J2 (5) 

m>n \ r m ~ r n\ 

Here qj is a charge of a j-th fragment. The forces Fj for SC are time-scaled as 

where pj n = pj — p n . Substituting in the Newton equation of motion one finds the 

following relation 

d 2 (j)(t) _ 1 ^ Pin , 

,. 2 m jPj = —77^ 1^ Wn — - (7) 

dt 2 0W 2 ^y Pin 

It is easy to see that it can be satisfied only if two conditions are fulfilled. Firstly, the scaling 
function should satisfy an equation 
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d 2 <P{t) = a_ , , 

dt 2 0(t) 2 ' 1 ' 

where a > is a time-independent constant which is discussed in detail below. One obviously 
recognizes in ([§]) the equation describing a one-dimensional motion of a particle with unite 
mass and unite charge in the attractive Coulomb field created by the charge a. Secondly, 
the validity of (|7|) needs that the vectors pj satisfy the following system of equations 

apj = Sj = -— Wk^r, (9) 

m j k^j Pjk 

They state that accelerations of each fragment ctj is proportional to its coordinate vector at 
the initial moment of time. Eqs.(^|) are shown to arise as conditions which are necessary for 
existence of SC. It is easy to see that they provide sufficient conditions as well. To verify this 
statement let us assume that we have a solution of (|9|). Then we can consider a trajectory 
with the following initial conditions. Firstly, we can choose initial coordinates as r}(£o) = Pj- 
Secondly, we can always choose initial velocities be proportional to coordinates 

^ =/Jr,«b), (10) 

where (3 is some positive constant which depends on the energy, (3 ~ \/E. From (^) we find 
that accelerations at the initial moment of time are also proportional to coordinates 

= aFifa). (11) 

Thus for the considered trajectory both the velocities and accelerations linearly depend on 
coordinates at the initial moment of time. Combining this fact with the Newton equations 
of motion we conclude that the velocities (and accelerations) remain to be proportional to 
the coordinates for any moment of time 

dfj{t) 



dt 



mm- (12) 



Here j3{t) is some positive function, (3{to) = (3. Integrating ([12]) we conclude that the time 
variation of distances does exhibit scaling condition (0), in which <p(t) = (3(t). 



This discussion shows that the SC (H) exists if and only if Eqs.(H) are satisfied. There 
are N vector variables pj, j = 1,2, ---N and one scalar variable a in these equations. 
Obviously not all of them are independent because there are seven transformations which 
do not change the given SC. Three of them correspond to shifts of the SC center of mass. 
Three others describe rotations of the SC as a whole. One more transformation describes 
the overall scaling of SC 

P3^^ = Xpj, j=l,2,---N, (13) 
a — > a' = X' 3 a, (14) 



with A > 0. According to Eq.(|8]) the scaling of a ($L4\) should be accompanied by a corre- 
sponding scaling of <f)(t), namely <ft(t) — > <j)'(t) = A~ 1 </)(t). Notice that the latter transforma- 
tion can be interpreted as a shift of the initial moment of time 

to ^4 (15) 
where according to Eq.(^|) t' Q should satisfy 

4t{Q = x-Wo) = 1- (16) 

It is easy to see that Eqs.(f|) remain invariant under the discussed above seven transforma- 
tions, i.e. the shifts, rotations and scaling, allowing one to consider them as a set of 3iV — 7 
equations for 3iV — 7 independent variables. When solving these equations it is convenient to 
treat a as a constant parameter which governs the overall scale and can be chosen arbitrary 
(for example a = 1). 

At SC the system Hamiltonian 

N v 2 df- 



is reduced to 



, 2m j ' I J J (II 



where 



N 



(19) 




(20) 



Clearly the Hamiltonian fll8|) describes the one-dimensional motion of a particle with the 
mass M. and unit charge in the attractive field of Coulomb center with the charge — Qq. 
The corresponding equation of motion is given by the considered previously Eq. @ in which 
the constant a proves to be equal to 

« = & < 21 > 

The interesting for us physical events take place if there is enough Coulomb attraction in 
the system. That is why we suppose that the effective Coulomb charge Q is attractive, 
Qo > 0, resulting in positive value of a. 

Eqs. (|T9D , (|20"D , (|2T|) show that arbitrary scaling of a can be compensated for by the cor- 
responding scaling of coordinates pj. This fact agrees with Eqs.([l3|),(|l4]). 

The scaling function <fi(t) is defined by straightforward integration of (H) 



where E is the system energy. Combined with the initial condition 0(to) = 1 this fixes the 
scaling function unambiguously. 



independent of the chosen coordinate frame. To see this more clearly let us introduce grand 
vectors in the 3N dimensional configuration space. The grand vector r(t) = (fi(t), ■ ■ -rjv(£)) 
defines the time-dependent coordinates, the vector p = (pi, • • ■ Pn) gives the initial coor- 
dinates and a = (ai,---ajv) is the vector of accelerations at the initial moment of time. 
We employ bold type to distinguish such a vector from the conventional vector in space. 
Eqs. ([!]), © allow the following presentation 




(22) 



It is important to emphasize that Eqs. ([[]), (§) present the idea of SC in an invariant form 
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r(t) = <f>(t)p, (23) 
ap = a. (24) 

Obviously, these relations between 3iV-vectors do not depend on a reference frame. This 
shows that the scaling coordinate (f)(t) is described in an invariant way. 

It has been presumed by previous authors that some coordinate which describes frag- 
mentation is to be singled out and the potential extremum point is to be found for the fixed 
value of this "break-up coordinate" . The latter has been chosen in most cases as the system 
hyperradius |3§ @ [@ @ defined as R 2 = E?= 

1 rriir i . In the hyper-coordinate reference 

frame the potential energy 

V=<M (25) 

is proportional to hypercharge C(u) which depends on a set of hyperangles u = 
(ui, • • - U3N-7)- It is eas Y to verify that definition of the SC (|4j) in the hyperspherical 
coordinates is reduced to 

= °- (26) 

which shows that a SC is a saddle-point of the hypercharge C(uo). The function <p(t) in 
hyperspherical coordinates is proportional to the hyperradius (j)(t) = R/Rq, where Rq is 
the initial value of the hyperradius. The effective charge Qq and the effective mass M. can 
be expressed in terms of i? an d the hypercharge M. = R%, Q = —C /Ro where C is 
the hypercharge evaluated for SC. Description of the system in the hyper-coordinates has 
a long tradition and list of achievements, see for the example recent calculations of the 
three-electron atom in hyperspherical coordinates [fiO~| |. However, generally speaking, these 



coordinates do not possess fundamental advantages over other coordinate frames for the 
fragmentation problem. 

Another well known reference frame provide Jacoby coordinates used in the approach 
developed by Feagin [51]. For the multiparticle fragmentation the choice of the "break-up 



coordinate" is not obvious and some special procedure was developed for its construction 
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29| f3~0|j |L9] . It becomes the more sophisticated problem the more complicated the system 



is. 



In conclusion of this Section it should be stressed once more that our approach provides 



an invariant definition for the idea of SC which is given in (23), (E 



III. SMALL DEVIATIONS FROM SCALING CONFIGURATION 



Assuming that the function <p(t) is defined as described in the preceding section, we 
switch from fj to the new coordinates Sfj 



fj = <j>(t) pj + 5fj 



(27) 



which have an obvious meaning of deviations from the SC. Presuming that these deviations 
are small, we write down linearized classical Newtonian equations for Sfj(t) as 



rrii 



d 2 Sri 



N 



d 2 



E v a 5r < 
=1 



<t>(t)>W ' 

E 



d Pi d Pj "n \Pm - Pi 

These equations of motion are generated by the time-dependent Hamiltonian function 

1 



(28) 
(29) 



N 



1 N 5v 2 

2 ^ mi W) 3 



E VtjSfi-Srj, 



St ' 

8% = m 3 ^. 



(30) 



It is convenient to introduce scaled deviations and related momenta 7?,- as 



J- 



(31) 



since this allows us to factor out the time dependence in the Hamiltonian: 

N ^2 i o n N 



5H 



(32) 



The derivation of this formula could be traced via a quantum mechanical analogue of the 
problem (which for many readers nowadays is more convenient than the pure classical consid- 
eration). In quantum mechanics the transformation rules for momenta and the Hamiltonian 
follow respectively from the formulae for the partial derivatives 
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k ~ m ~ V4 di; { d dt) sfj ~ U) 6 + 40(t)V4 * ■ 4- ( 33 ) 

In (^) we use a symmetrized representation which should be employed in the quantum 
version of the formulae (the latter also implies a corresponding gauge transformation for the 
wave function). 

From ( p2"D one obtains 



<ty 2(E<j>(t) + Qo) 

dt=i — m — (34) 

which becomes time- independent for E = 0. In this case the time-dependence is exactly 
factored out in the Hamiltonian ([32]) justifying the choice of the scaling fl3T|). This implies 
that the original non- stationary problem becomes stationary provided one replaces time t 
by an effective time r. A relation between t and r in differential form is 

dr = (j)(t)-*l 2 dt. (35) 

For some applications it is necessary to keep the energy dependence of the trajectory. For 



these cases a convenient technique has been developed recently by Kuchiev [17]]. We have 
applied it to the case considered and verified that it results in the same threshold indexes 
as the ones obtained below by the stationary approach. 

The Hamiltonian describing propagation in the effective time (53) reads 

1 N 7? 2 1 a N 

SH r = 2 E ± + o 5>«S§ + 2 £ (S • *i + *i -5) > (36) 

j=l i ij j=l 

- = -|V5f- (37) 

The Hamiltonian 5if T ( p6|) is quadratic in coordinates and momenta thus describing a 
set of harmonic oscillators or inverted oscillators. This shows that our goal is to describe 
the behavior of the system in terms of these oscillators and inverted oscillators. Before 
proceeding we modify our notation. The set of components of the displacements vectors Sfj 
(j = 1,2,... N) comprise 3A^-dimensional grand vector St. In this formulation, for instance, 
Vij corresponds to grand 3iV x 3iV square matrix denoted below as V. We introduce also 
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3iV x 3N unit matrix I and the diagonal matrix K of the same size with diagonal elements 
corresponding to inverse mass l/rrij of each particle. 

This notation takes into account an obvious fact that the total number of all modes 
coincides with the number of degrees of freedom in the system (k = 1, 2, ... 3N). There 
are however seven particular degrees of freedom: translations, rotations and the scaling 
transformation. They do not change the shape of a SC and do not describe a deviation from 
a SC. These degrees of freedom may be called the collective modes. They obviously should be 
considered separately from the oscillating modes which describe deviations from the SC. In 
order to distinguish the collective modes one can use the following interesting property. All 
collective degrees of freedom are described by the eigenvectors of the grand matrix KV with 
particular eigenvalues. Firstly, the three modes which correspond to the system translations 
in space have obviously zero eigenvalues. Secondly, the modes corresponding to rotations of 
the system in space have eigenvalues equal to Qo/JA, as shown in Appendix. There are three 
such modes in general case, while for a linear SC there are only two modes. Thirdly, the 
mode corresponding to the scaling transformation Eqs.([T^),(I3) has an eigenvalue — 2Q /Ai, 



as also shown in Appendix. Using these eigenvalues one can separate the collective modes 
either from the very beginning, or at the end of calculations. 

There is another useful for applications way to separate the collective modes. For trans- 
lations and rotations the separation can be fulfilled by conventional methods choosing ap- 
propriately the coordinates, as is demonstrated in a number of examples below. Separation 
of the scaling mode can be achieved with the help of the operator of projection on this 
mode P and the complementary projection operator Q = I — P. The operator P is readily 
constructed from the unit vectors fij = Pj/pj which define the shape of SC: 

Pij = Hi ■ fij. (38) 

Thus all seven collective modes can be easily identifies and separated using any of the two 
techniques described above. 

Some modes in the 3N — 7 subspace orthogonal to the collective modes are stable and 
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describe small oscillations around SC; the related oscillating frequencies Uk are real. The 
object of our major interest is unstable modes with imaginary oscillating frequencies. It is 
shown below that unstable modes exist for any SC. It is convenient to introduce for unstable 
modes a parameter au = iuik (Rea^ > 0). In order to find the oscillating frequencies one 
can presume a harmonic time-dependence of the coordinates £ and momenta n 

£ = exp(iut) E, 7r = exp(iut) II, (39) 

where 3 and II are time-independent grand vectors. The Hamiltonian equations of motion 
give 

iw3 = Kn + a3, «n = -V3-an, (40) 

where a (|37| ) is a scalar coefficient. The latter equation could be written also as 



1 

rii 


( a K\ 


[I] 


H 




H 



(41) 



Excluding the grand vector II one comes to the eigenvalue problem for the square of fre- 
quency u) 2 

(w 2 + a 2 )3 = KVS, (42) 

or, in the symmetrized form 

(a; 2 + a 2 )H = K 1 / 2 VK 1 / 2 H, K 1 / 2 S = S . (43) 

Denoting a set of eigenvalues of the matrix KVas^, h — 1, 2, • • • , 3N, we obtain 

u 2 k = v k - a 2 , (44) 
«fc = \Ja 2 ~ v k . (45) 

This formula shows how the oscillation frequencies depend on the eigenvalues of the matrix 
KV. 
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Let us verify now that a SC is always unstable. With this purpose let us show that the 
matrix KV always possesses negative eigenvalues which describe instability. Consider the 
trace of the grand matrix V 

TrV = £-^— £-^-=£A^. £-^- = 0, (46) 

■ UpjUpj m>n \p m p n \ ■ m>n \Pm Pn\ 

which vanishes since the Coulomb potential satisfies the Laplace equation 

1 _ t _ 

A ft \n -n \ = Pj^Pn- 
\Pj Pn\ 

It is easy to see also that ( f46"D results in Tr (K V) = which means that 

Tr(KV)=5>k = 0. (47) 

k 

We see that the spectrum of the matrix KV always contains both positive and negative 
eigenvalues. This fact in itself is not sufficient to make a statement about instability because 
the trace ( fT7|) includes contribution from collective modes which do not change a shape of 
the SC. However, it is easy to exclude collective modes. Remember that the eigenvalues 
corresponding to translations are zero, rotations give eigenvalues Q/A4, while the scaling 
transformation provides the eigenvalue — 2Q/Ai, see Appendix. The sum of eigenvalues of 
collective modes is 

!QI M. in general case 
(48) 
for linear SC. 

Subtracting this result from ([|7]) we find the trace of the matrix K V in the subspace 
orthogonal to the collective modes 

{— Q/M. in general case 
(49) 
for linear SC. 

Since this trace is non-positive, we conclude that the matrix K V inevitably possesses nega- 
tive eigenvalues which describe deviations from the SC. This shows that any SC is unstable. 
This property is closely related to the fact that harmonic functions, i.e. those which satisfy 
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the Laplace equation, cannot have maxima or minima. Remember also the Earnshow the- 
orem well known in electrostatics: stable equilibrium is impossible for systems where the 
Coulomb forces are operative. Although SCs describe expanding non-static configurations, 
conclusion about inevitable instability remains valid in this case as well. This fact can be 
interpreted as a dynamic analogue of the Earnshow theorem. 

IV. QUANTIZATION OF DEVIATIONS FROM SCALING CONFIGURATION 

AND THRESHOLD INDICES 

Previous section reduces description of small deviations from SC to the set of coupled 
harmonic oscillators which could be quantized straightforwardly. This procedure provides 
the 'energy' levels 



Here the first subscript k — 1, 2, • • • , 3iV — 7 indicates the mode, and = 0, 1, • • ■ shows a 
number of quanta in this mode. For a given set of the quantum numbers {n^} the system 
wave function is given by 



where we omit the common time-dependent phase factor. The wave function is prepared at 
some initial moment to by preceding strong interaction of all fragments. In the Wannier- 
type approach it is presumed that these processes depend smoothly on the energy E. Hence 
they do not influence the form of threshold law and thus could be effectively excluded from 
consideration; it is sufficient to consider only t > to domain. 

For unstable modes the 'energies' e^ nk are complex-valued which leads to the loss of 
probability in the expanding SC. This should be interpreted [[34]] [|35| as sliding from the 
potential saddle in multidimensional configuration space that eventually leads to formation 
of bound states of two (or more) fragments. Such an outcome implies that the related part 




(50) 




V I Jt J \ l Jto <p\tp 



(51) 
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of probability is lost for the process of complete system fragmentation which is an object of 
our study. Cross section of the latter is proportional to the survival probability 



P, 



{ n k} 



{n fc }|i— >oo 



2 



exp (-V2M T r 7 ^ ^ d<p(n k + -)) , (52) 



where summation over k runs over all unstable modes. Note that the original quantum 
problem is stationary. The time t in ( |52"D plays a role of an effective variable which describes 
scaling of the system in accordance with (0). Small deviations from SC are described 
quantum mechanically. Our treatment generalizes to multimode case the scheme developed 
by Kazansky and Ostrovsky Q for the two-electron escape (see also Ref. ||12|| ; some of ideas 



used were elaborated also by Watanabe 0). Note that the cited paper ||] provides also a 
description of deviations from pure power threshold law, but we do not pursue this point 
here. 

For our objectives it is sufficient to note that fl5"2"|) has a form of a product of contributions 
coming from each individual mode, hence the threshold law of interest is 

a ~ P {nk} ~ (53) 

fi = J2fi knk . (54) 

k 

The partial threshold indices fi krik stem from the 'eigenfrequencies' of unstable modes being 
related to the negative eigenvalues v k < of the K V matrix 



fik 



M 9 1 

v k + 777 7 



n k + -) . (55) 



2Q 16 4_ 

Small positive values of v k formally could also lead to real fi knk , but have to be discarded. 
Obviously, if some (imaginary) 'eigenfrequencies' are iVfc-fold degenerate, the related contri- 
butions appear N k times in the sum (0). In principle the wave function is a superposition 
of terms corresponding to various sets of quantum numbers {n k }, since all of them are pop- 
ulated by the processes in the inner interaction domain. Clearly, the threshold law is defined 



by the least possible values of n k E| which are equal to zero unless the symmetry consid- 



erations forbid this choice, as exemplified in the next paragraph. If the initial SC is scaled 
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by the factor A, see Eqs.flT3|), (0), then KV~ A" 3 , M ~ A 2 , Q ~ A -1 , but the threshold 
indices Pkn k , as anticipated, remain scale-independent. Note also that the threshold index 
is invariant under simultaneous scaling of all charges or all masses in the system. 

In the original Wannier problem two electrons escape from infinitely heavy atomic core 
with the charge Z. The configuration found by Wannier gives the simplest example of SC 
in which the electrons reside at equal distances p and in opposite directions from the core. 
The motion is unstable with respect to the stretching mode which is separated from the 
(stable) bending mode. Thus it is sufficient for our purposes to consider motion of electrons 
along the line passing through the core. This motion is described by two coordinates and 
the matrix V takes the form 

( ~2Z ' 1 1 ^ 



p 3 



V 



-2Z 



4/ 



and <2o = \2Z — |J p -1 , M. = 2p 2 (we use atomic system of units, K = I). The eigenvalues 
of V are V\ = —2Z/p 3 and v 2 = (^—2Z + |) /p 3 . The eigenvalue v 2 is seen to coincide with 
— 2Q/A4. Hence it corresponds to SC expansion (see Appendix) and should be discarded. 
The eigenvalue v± upon substitution into (|55|) reproduces the well known result 



'100Z - 9 
AZ- 1 



1 

Til H — 
1 2 



(56) 



The choice n\ — provides the famous Wannier law valid for 1 S symmetry of the final 
two-electron continuum state, whereas ri\ = 1 corresponds to the threshold law for 3 S e (and 
3 P e ) symmetry 

Feagin and Filipczyk and Poelstra el al [ 30 1 put forward another formula for the 
threshold index in the multimode case. According to it the Wannier index is (N — 2) times 
larger than (|54l). The factor (N — 2) is described as a "phase space factor for (N — 1) 
outgoing particles" being justified by the reference to the earlier paper by Feagin We 
were unable to find the derivation of such a factor in the cited paper; anyway it deals only 
with the conventional N = 3 case where the factor (N — 2) is insignificant. Our treatment 
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provides purely dynamic approximation for the wave function and does not leave any room 
for the statistical arguments. The other aspects of relation between dynamic and statistical 
threshold laws are discussed in Section VI. 



V. PARTICULAR SYSTEMS 

In practical applications of our scheme the less obvious part corresponds to finding SCs. 
Numerical solution of the set of non-linear equations @ could be cumbersome and implies 
reasonable initial guess. The question whether all the solutions are found is even more 
difficult. In reality one has to appeal to intuitive reasoning and to limit search to some 
symmetrical configuration. This allows one to effectively reduce the number of equations 
(|) to be considered. Since the initial step of finding SC in most cases could not be done 
in closed form, we do not pursue the goal of obtaining analytical formulae, but resort to 
numerical calculations which are performed using the Mathematica program. We find it 
easier to avoid preliminary separation of rotational and translational coordinates, since they 
could be easily distinguished in the eigensystem of the complete matrix K V. Moreover, the 
known eigenvalues of this matrix corresponding to rotations (see Appendix) provide a good 
test for consistency of calculations. 

The systems practically accessible nowadays in atomic physics are not very diverse, 
consisting of several electrons and positrons in the field of heavy (positively charged) atomic 
core. Since three-particle systems (such as A +z + 2e~ or A +z + e~ + e + ) are already studied 
in great detail [|J ]38| j|5| [[4(J j|7| (see also references in the Introduction), we start 
from the four-particle systems. We do not impose any symmetry constraints on the system 
state thus presuming that = for all modes contributing /i (jj4]). 

A. Three-electron escape from the charged core 

The system A +z + 3e was thoroughly investigated by Klar and Schlecht and Grujic 



2J|. They considered a configuration of electrons forming an equilateral triangle with in- 
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finitely massive core in the center, which is obviously a SC. The out-of-plane motion is 
separated. It corresponds to stable modes and does not affect the threshold law. The in- 
plane motion is described by six coordinates of electrons, or by four 'oscillatory' modes plus 
uniform expansion of SC and its rotation. The eigenfrequencies obtained by us, as well as in 
the cited papers, are pairwise degenerate due to SC symmetry. One pair corresponds to sta- 
ble motion and the other pair to unstable motion. The latter pair produces two equal terms 



in the sum (|54|) . Klar and Schlecht [22[ and Grujic [24| succeeded in deriving analytical 



expressions for the Wannier index |48|]. In this paper we do not pursue analytical formu- 
lations but check that our numerical results coincide with those cited by Grujic, namely 
H = 2.82624 for Z = 1, // = 2.27043 for Z = 2, // = 2.16196 for Z = 3, etc. The experiment 
for electron impact double ionization of atoms [Z = 3) seem to agree with the threshold law 



The two pairs of modes discussed above are already well known. Combined with rotation 
and scaling expansion they represent a complete set of six in-plane coordinates. Since the 
number of modes is a physical parameter which is independent on the theoretical technique 
used, we do not see any possibility to obtain some additional unstable modes which would 
lead to another Wannier index and thus to the complementary threshold law as announced 
by Feagin and Filipczyk |29| (in fact our conclusion could be drawn from the paper by Grujic 



who used the plain Cartesian coordinates whereas less transparent treatment by Klar 
and Schlecht |22|] is based on hyperspherical coordinates). Since no details of analysis by 
Feagin and Filipczyk |29| were ever published, more detailed discussion of this issue is not 



possible. 

B. 2e~ + e + escape from the charged core 

The plausible symmetric SCs for the system A +z + 2e~ + e + were considered by Poel- 
stra et al [[HJ (note that the calculations in this paper were carried out only for Z = 1). 
They comprise two different linear arrangements and one plain configuration [|49|| . All these 
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configurations belong to SC and therefore can be easily handled by the developed above 
technique. We consider below these SCs successively 



1. Linear configuration L a 

Let the frame origin is placed into infinitely massive core having the charge Z. The 
coordinates of two electrons and positron are respectively X\, x 2 , x 3 ; all of them are positive. 
It is convenient to introduce two dimensionless parameters x = ri/r 3 and y = r 2 /r 3 (0 < 
x < 1 < y) which have to satisfy the system of equations obtained from @ 

m 3 Z/x 2 - 1/(1 - x) 2 + 1/(3/ - x) 2 



mi -Z + 1/(1 -x) 2 - 1/(1 -y) 2 
m 3 Z/y 2 + 1/(1 -yf-l/{y- X f 



x. 



y (57) 



m 2 —Z + 1/(1 — x) 2 — 1/(1 — y) 2 

(equations are presented for more general case when all light particles have different masses 
rrii, while the core remains infinitely heavy). 

2. Linear configuration L^ 

A distinction from the previous case is that the coordinate of one of the electrons is 
negative (x% < 0). The system of equations defining SC is somewhat different (y < < x < 
1): 

m 3 Z/x 2 - 1/(1 - xf - l/(y - x) 2 _ 

x, 



mi -Z + 1/(1 -x) 2 + 1/(1 -y) 2 
m 3 -Z/y 2 - 1/(1 -yf+l/(y- X ) 2 



y. (58) 



m 2 —Z + 1/(1 — x) 2 + 1/(1 — y) 2 

For both linear configurations the bending modes are stable. There are two stretching modes 
for each configuration, both being unstable. The results of our calculations are summarized 
in table 1. For Z = 1 parameters x, y and partial threshold indices \i\ and fi2 coincide with 



those obtained by Poelstra et al j30[; our threshold indices /i are less by a factor of 2, as 



discussed at the end of Section [IV|. Notice a non-trivial behaviour of the parameters with Z: 
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for instance, in SC L a , x and /xi increase with Z, whereas y and \i2 decrease. The threshold 
index \i increases with Z which is opposite to the well known behaviour for the simplest 
system (A z + 2e) and for 3e escape where /i diminishes as Z grows (see more discussion in 
the Section |V]). 



3. Plane configuration P 

The symmetric plane configuration is conveniently characterized by two angles a and 
(3: a is an angle between two lines which join the ion with the positron and one of the 
electrons, while (3 is an angle between two lines which join the positron with the ion and 
with one of the electrons. In the considered plane configuration the two electrons are located 
symmetrically, which means that their locations mirrow each other under reflection in an 
axes which joins the ion and the positron. This makes the angles a and (3 be identical for 
both electrons. From (|9|) we deduce the system of equations 

m + sin 3 7 i^Z sin 2 a — - sin a + sin 2 (3 cos 7^ = 
= m_ sin 3 (3 (2 cos (3 sin 2 7 — Z sin 2 an , 
sin 2 (3 sin 7 = - cos a (7 = it — a — j3), (59) 

where the masses of light particles with negative (m_) and positive (m + ) charges generally 
could be different. The results of calculations are presented in table 1. For Z = 1 the angles 



a and (3 coincide with these extracted from the paper by Poelstra et al pOfl . However, 
the difference between the threshold indices is drastic. Poelstra et al had found a single 
unstable mode which corresponds to our partial Wannier index /i 2 . Our calculations give 
two unstable modes, similarly to the case of 3e escape (in the latter case the modes were 
degenerate due to a symmetry which is absent for the system under consideration). The 
reason of this disagreement remains unclear. 

The plane SC P governs the threshold behaviour, although it provides the threshold 
index /j only slightly less than the linear configuration L^. 
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C. Four-electron escape from the charged core 

Basing on the symmetry considerations we analyze three configurations: linear, plane 
and 3D SC. It could be shown rigorously that for symmetric linear arrangement SC does 
not exist for all values of Z, i.e. Eqs.(Q) have no solution. 

1. Plane configuration P 

In the plane configuration the electrons are located in the apexes of a square; the core 
lies in its center. The out-of-plane motion is separated and corresponds to stable modes. 
For in-plane motion in general case we find one non-degenerate and one doubly-degenerate 
unstable modes (table 2). For the particular case Z = 1 an additional non-degenerate mode 
becomes unstable. 

2. 3D configuration V 

SC describes the electrons located at the apexes of tetrahedron. We find a single triply- 
degenerate mode (table 2). Interestingly, the threshold index fi proves to be quite close for 
plane and 3D configurations, although 3D SC provides somewhat lower value of \x and thus 
governs the threshold behaviour. As Z increases, the relative importance of electron-electron 
interaction decreases and \x approaches the value \x = 3 which corresponds to non-interacting 
electrons. 

The smallest practically attainable value of the charge seems to be Z = 2. It could be 
realized via triple ionization of negative ion by electron impact. However, theoretically the 
case Z = 1 proves be very interesting due to unusual properties. In this case the threshold 
index becomes much larger than in other cases, particularly for the plane SC. This is due to 
a small value of the 'charge' Qo in this case. Another interesting feature is an appearance of 
an additional unstable mode in the plane SC. An analysis of the eigenvector E shows that 
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it corresponds to the out-of-plane motion. Namely, a pair of electrons lying on a diagonal 
of the square shifts upwards, whereas another pair shifts downwards. 

The tetrahedric configuration was considered earlier by Grujic |26[ who obtained ap- 
proximate analytical expressions for the threshold indexes. The partial threshold indexes 
obtained by him reveals only an approximate degeneracy. The numerical results for \i are 
in reasonable agreement with our data. 



D. 3e + e + escape from the charged core 
We failed to find a symmetrical plane SC for this system. 



1. Linear configuration L 

Linear SC corresponds to alternating positive and negative charges. Let us locate the 
origin at the heavy ion and call by x\ > 0, £2 < 0, £3 > x\ > locations of three electrons, 
and by £4, X\ < £4 < £3 location of the positron. Then the considered configuration can be 
characterized by three parameters: x — X\/x^ y = £2/^4, z = £3/0:4 ((y < < x < 1 < z). 
They have to satisfy a set of equations which follow from (fj|) 



1 



1 



1 



£ 2 (x — y) 2 (1 — x) 2 (z — xf 
Z 1 1 1 



£ 



y 2 (x-y) 2 (l-y) 2 {z-yf 
Z 1 1 1 

^7 + 7 7T7 + 



(z — x) 2 (z — y) 2 {z — 1)' 



z 
z 
z 



1 



1 



1 



(1 


-x) 2 (z-l) 2 


(i-y) 2 . 




1 1 


1 


(1 


-x) 2 ' (z-l) 2 


(1 -2/) 2 J 




1 1 


1 



[l-xf (z-l) 2 (l-y) 



(60) 



The parameters of SC and the Wannier indices are shown in table 3. 



2. 3D configuration V 

The symmetrical 3D configuration arises when the three electrons form the equilateral 
triangle while the ion and the positron are located up and down the plane of the triangle 
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on the perpendicular to the triangle plane which crosses its center. This configuration is 
characterized by two angles a and (3 defined similar to the case considered in Section B3. 
Namely, a is an angle between two lines which join the ion with the positron and with one 
of the electrons, while /3 is an angle between two lines which join the positron with the ion 
and with one of the electrons. These angle are difined by equations similar to 

m + sin 3 7 ( Z sin 2 a — —j= sin a + sin 2 (3 cos 7 J = 



m_ sin 3 (3 (3 cos (3 sin 2 7 — Z sin 2 a 



sin 2 (3 sin 7 = —= cos a (7 = 7r — a — f3). (61) 

v3 

We have found two doubly-degenerate and one non-degenerate unstable mode as shown in 
table 3. The threshold law is governed by 3D SC V. Note that the threshold index grows 
with Z. 

E. Five-electron escape from the charged core 

1. Plane configuration P 

In the plane SC the electrons are located in the apexes of a equilateral pentagon; the 
core lies in the same plane. In the in-plane motion we have found two doubly degenerate 
unstable modes (table 4). For Z = 2 an additional pair of unstable modes appears. 

2. 3D configuration V 

Here three electrons lie in the apexes of equilateral triangle with the core in its center. 
On the perpendicular to this plane, above the plane and below it, another pair of electrons is 
located symmetrically. The SC can be characterized by the angle a between the line which 
joins out-of-plane electron with the core and the line which joins it with in-plane electron. 
The angle is defined by the equation 
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—j= + 2 sin 3 a — Z = tan a ^3 sin 2 a cos a + - tan 2 a — Z tan 2 o^j . (62) 

Quite unexpectedly, a proves to be very close to 45°, exhibiting weak dependence on the 
core charge Z (table 4). This means that in-plane and out-of-plane electrons are located at 
almost the same distance from the core. The 3D SC generates somewhat lower values of \i 
than the plane SC thus governing the threshold behaviour. However, the difference is quite 
small. This feature is common to that found above for the four-electron case. 

3D configuration for five-electron system was considered previously by Dmitrievic et al 



27| . However, the equation derived for SC angle a differs from (p2|). 



F. Fragmentation in two pairs of identical particles with opposite charges 

In this subsection we consider fragmentation into the final state 2 X^ z + 2e, where X^ z 
is a positively charged particle with charge Z and mass m (all results below hold if the 
electrons are replaced by any other charged particles; then Z and m have the meaning of 
ratio of charges and masses respectively). In the applications considered above the zero 
eigenvalues of the matrix V do not emerge due to the presence of infinitely massive core. 
In the (2X+ Z + 2e) system such modes are present. Another distinction is that for equal 
masses of leptons in previous applications we have always had K = I and v k have been the 
eigenvalues of the V matrix. Now we have to diagonalize the complete matrix K V. Both 
these features do not create substantial difficulties. 

From the symmetry considerations it is clear that a shape of the SC is a rhombus with 
the angle 2a at the apexes where the particles are situated. The single SC parameter a 
is defined by the equation 

8Z ^— = m ( 8Z - -\-] ■ (63) 



cos 3 a V sin 3 a 
which follows from (|9|). Several examples are shown in table 5. The simplest practical 
realization is the complete fragmentation of H2 molecule by photons where a is close to 
30° in agreement with Feagin and Filipczyk [^] and the threshold index proves to be huge. 
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Apparently this threshold behaviour could not be observed in experiments |pOfl . Another 
feasible realization with moderate Wannier index is ionization of negative positronium ion 
by positron impact (Z = 1, m =1). 



We fail to find the linear configuration discussed by Stevens and Feagin [31 



VI. DISCUSSION AND CONCLUSION 

This paper formulates the idea of the SC. Defined by (p]), the SC is shown to arise 
when a nonlinear set of Eqs.(BI) is satisfied. Propagation of the system in the vicinity of 



SC configuration governs the threshold law which is found in Eqs.(|53D,(p4|) and (pq) . These 
results permit direct practical calculations of the threshold index \i for any system. 

In many cases the threshold laws in quantum mechanics can be deduced from general 
considerations without dynamical treatment. For instance, the break-up cross section with 
N fragments in the final state and a short range interaction between them could be estimated 
from simple phase-space volume (i.e. statistical) arguments as 



a s „ £§("-i)-i. (64) 

If one presumes that all fragments ("electrons") are attracted by Coulomb forces to one 
fragment ("core"), but the interaction between the "electrons" is negligible, then the phase 
space arguments could be easily modified to give 

a c ~ E N ~ 2 . (65) 

In case of repulsive Coulomb interaction with the "core" (but still without other inter- 
fragment interactions) the cross section at the threshold becomes exponentially small, as 



obtained, for example, by Geltman |HJ in his calculations for atom ionization by positron 
impact with all correlation neglected. The threshold behaviour changes to ~ E 3 ^ 2 [^] if 
one employs the so called 3C wave functions for the final continuum state. However, these 
functions do not ensure proper description in the near-threshold domain. 
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If one aims to obtain a correct threshold law for the Coulomb system, then the interaction 
between the fragments, i.e. the particle correlation is to be taken into account. This makes 
the phase-space arguments insufficient, but requires dynamical treatment as it was originally 
done by Wannier for the simplest system. In this paper we employ the most simple 
theoretical apparatus presenting the essential equations in an arbitrary coordinate frame. 
They remain valid, in particular, in the simplest single-particle Cartesian coordinates. 

As discussed in Section [TV], Poelstra et al [f3C§ suggested another formula for the Wannier 
index which differs from our Eq. fl54j) by the extra "phase factor" (N — 2). This discrepancy 
remains hidden when one restricts consideration to the case of two, three, or four electrons 
receding from the positively charged core. In these cases the unstable mode proves to be 
respectively non-degenerate, doubly- and triply-degenerate. Thus the degree of degeneracy 
in these cases coincides with (N — 2). This fortuitously allows one to replace the summation 
over degenerate modes implied by formula (|54] ) by multiplication over the factor (N — 
2) which corresponds to the formula by Poelstra et al |[30| . However this coincidence is 
accidental and misleading. It is broken, for instance, by variation of charges and masses of the 
constituent particles which violates SCs symmetry and hence lifts the modes degeneracy, or 
by considering larger numbers of particles N (simply because possible degrees of degeneracy 
are restricted by properties of the point groups in 3D space). For five electrons receding 
from a charged core only doubly degenerate unstable modes were found above. 

Physically it is clear that if the charge of the core Z in the system A z + (N— 1) e becomes 
bigger, then the interelectron correlations should become less important and the threshold 
law should approach the value obtained from the phase-space arguments, i.e. — > (N — 2) 
as Z — > oo. This conclusion is supported by all examples considered. Moreover in all these 
examples one can note that: 

• the number of unstable modes accounting for their degeneracy (i.e. the number of 
terms in the sum (|5lD) is equal to (N — 2); 

• each partial Wannier index fi^o (11) tends to unity from above as Z increases. 
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Apparent exception from the first rule is an emergence of an additional unstable mode 
in the plane A +z + 4e SC for Z — 1. However, this SC provides [i larger than 3D SC and 
therefore it does not govern the threshold behaviour. Note that although these properties 
are physically very natural, it is not clear if they can be proven rigorously from the first 
principles. An additional observation is that the electrons in SC tend to be distributed 
uniformly on the sphere, even when the corresponding perfectly symmetrical body does not 
exist (see five-electron case above). For large number of electrons in the field of the core 
several competing SC are found to produce very close threshold indices. Still, in all the cases 
considered the leading SC is found to be the three-dimensional one. 

These results hopefully should hold if the electrons are replaced by other (possibly dif- 
ferent) negatively charged particles. However, the situation changes drastically if one of the 
"electrons" is replaced by a particle of positive charge, for example, positron. It is essential 
that an additional repulsive Coulomb interaction appears in the system. If correlations are 
neglected then the cross section decreases exponentially as E approaches threshold. One 
could expect that although the true threshold law retains a power character for all values of 
Z, it tends to mock the exponential behaviour by increasing of /i value [53|. This property 



holds for all positron-containing systems considered above. The threshold index increases 
with Z quite slowly. In order to illustrate the later point quantitatively we cite results for 
A +z + 2e~ + e + system with very large values of Z (cf. Sec. [VjB): jj, = 9.4 for Z = 50 
(a = 5.60°, $ = 37.2°); // = 11.6 for Z = 100 (a = 3.95°, (3 = 37.8°). In general terms 
one can argue that a similar situation should arise when a system contains two or more 
positively charged particles and two or more particles with negative charge. Note that the 
properties of the partial Wannier indices /i^o are less straightforward: some of them could 
be less than unity and vary with Z non-monotonically. 

Large values of threshold indexes /i are unfavorable for an experimental observation of 
the threshold behaviour: close to the threshold the cross section proves to be too small 
to be observable, and for higher excess energies the intrinsic deviations from the threshold 
law become essential. An analysis of the energy domain where the threshold law holds is 
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beyond the scope of this paper. Still, we can note that for the electron-impact ionization of 
atoms or for double photoionization this domain is limited to few eV above threshold (for 



quantitative treatment within the Wannier mechanism see Refs. [ffl |T3]). For the positron- 



impact ionization the applicability domain is even less [O E| . As argued by Ihra et al |L8 



an agreement with experimental data could be substantially improved if the interaction of 
different modes in the deviation from SC is taken into account. Possibly some procedure 
to assess for the mode interaction could be developed also for the multifragment system; 
the present development provides a necessary first step for more advanced approaches. One 
could note also that even very large threshold indices could (quite unexpectedly) be useful 
for constructing formulae of interpolation character as shown in the recent paper by Rost 
and Pattard pO . 
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APPENDIX A: EIGENVALUES OF V MATRIX CORRESPONDING TO 
ROTATIONS AND TRANSLATIONS IN TIME 

If the V-particle system is rotated as a whole over infinitisemal (time-independent) angle 
S(p around the axis u, then the particle coordinates receive increments 

Sij v) = (0 x fj) Sip. (Al) 

The form of Newtonian equations of motion 
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d 2 f- BTI 



remains invariant under rotations. This implies that 

d 2 5r< JL 8 2 U 



™™ A = -Et^-^ (A 3 ) 



3 cfa 2 ^ <9r*j drj 

where 5fj = Srj. For SC one can use Eqs. (|A1|) and (P to get 



dt 2 dt 2 J 
Bearing in mind that according to ([18]) 

^ dt 2 M 
and using definition (|29D , we finally obtain 



°— -5r)'. (A4) 



(A5) 



1 AT n 

-E^^ W = ^^ W (A6) 
m 3 £i A4 

which means that the grand vector 5r ^ is an eigenvector of the grand matrix K V with the 
eigenvalue Qo/ ' M.. Generally there are three eigenvectors corresponding to this eigenvalue, 
but for a linear SCs only two independent rotations are possible. 

Consider now variation of the trajectory caused by shifting of time over an infinitesimal 
interval t — > t + 5t using similar technique. For the system in SC the particle coordinates 
are incremented in this case by 

5ff c) = vjdt = ^ pj St. (A7) 
etc 

The form of Newtonian equations of motion obviously remains invariant under the shift of 



the time variable. Therefore (|A3|) remains valid for 5fj = 5fj . An analogue of Eq. (|A4| 



.(SC) 

diiauic. j-iiciciuic ^|i~vo| y icinciiiio vaiiu iui uij — ui 

now reads 



dt 2 \dt ) dt 3 
Differentiating Eq . (|A5|) we obtain 
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dt 3 M 3 dt 
which finally brings us to 



rf 3 <#> 2So 1 # (A9) 



±gv jl «f°> = _aoi«f9 . (mo) 



Since the grand vector <5r ( - sc - ) is proportional to the grand vector p which defines the SC 
shape, we conclude that the latter vector is an eigenvector of the grand matrix K V with 
the eigenvalue — 2<2 /.M. 
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TABLES 



TABLE I. Parameters of scaling configurations and Wannier indices for A +z +2e +e + system. 



z 




SC parameters 




Mi 


M2 




SC L a 


1 


a; 


= 0.506100, 


y = 1.692952 




4.442178 


2.193945 


6.636123 


2 


X 


= 0.587468, 


y = 1.636629 




4.767141 


2.064237 


6.831377 


3 


X 


= 0.633155, 


y = 1.609587 




5.024502 


1.966884 


6.991386 


4 


X 


= 0.664214, 


y = 1.594313 




5.242328 


1.88483 


7.127158 


SC L b 


1 


X - 


= 0.441380, 


y = -0.677611 




2.577720 


1.025435 


3.603155 


2 


X - 


= 0.539724, 


y = -0.847969 




2.888492 


1.009213 


3.897705 


3 


X - 


= 0.594480, 


y = -0.949091 




3.193559 


1.005040 


4.198599 


4 


X - 


= 0.631720, 


y = -1.023071 




3.475766 


1.003244 


4.479010 


SC P 


1 


2a 


= 76.7338°, 


2(3 = 55.1969° 




1.884950 


1.562234 


3.447184 


2 


2a 


= 55.1741°, 


2(3 = 61.3793° 




2.045028 


1.793101 


3.838128 


3 


2a 


= 45.4233°, 


2(3 = 64.1787° 




2.206553 


1.972092 


4.178645 


4 


2a 


= 39.5138°, 


2(3 = 65.8916° 




2.351217 


2.123469 


4.474686 
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TABLE II. Wannier indices for A +z + 4e system. The numbers in parentheses indicate degree 
of unstable mode degeneracy. 



z 


Mi 


M2 


^3 


A* 


SC P 


1 


4.877419 


4.248225 (2) 


2.071837 


15.44571 


2 


1.356093 


1.273381 (2) 




3.902855 


3 


1.192808 


1.145660 (2) 




3.484128 


4 


1.132414 


1.099316 (2) 




3.331046 


5 


1.100871 (2) 


1.075346 




3.251563 


SC V 


1 


3.075960 (3) 






9.227870 


2 


1.257986 (3) 






3.773958 


3 


1.139795 (3) 






3.419384 


4 


1.095940 (3) 






3.287819 


5 


1.073040 (3) 






3.219120 
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TABLE III. Parameters of SCs and Wannier indices for A +z + 3e + e + system. The numbers 
in parentheses indicate degree of unstable mode degeneracy. 



z 




SC parameters 




Mi 


M2 




A* 


SC L 


1 


x = 


0.580448, y = - 


1.070391 


z = 1.627861 


4.41213 


2.30976 


1.06890 


7.79079 


2 


x = 


0.580448, y = - 


1.070391 


z = 1.627861 


4.73844 


2.13221 


1.03343 


7.90408 


3 


X = 


0.628772, y = - 


1.162883 


z = 1.602043 


4.99954 


2.01766 


1.02270 


8.03998 


4 


X = 


-- 0.661096, y = - 


-1.22438 


z = 1.587485 


5.21972 


1.92598 


1.01746 


8.16316 


SC V 


1 




a = 60.5698° , 


P = 32.2041° 


1.57584 (2) 


1.03194 


0.60493 (2) 


5.39348 


2 




a = 40.5400° , 


P = 41.7154° 


1.56354 (2) 


1.20043 


0.66302 (2) 


5.65356 


3 




a = 32.3675°, 


P = 44.9869° 


1.70957 (2) 


1.33711 


0.62771 (2) 


6.01166 


4 




a = 27.6668°, 


P = 46.7663° 


1.85129 (2) 


1.45098 


0.57327 (2) 


6.30011 
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TABLE IV. Parameters of SCs and Wannier indices for A +z + 5e system. The numbers in 
parentheses indicate degree of unstable mode degeneracy. 



Z SC parameters [i\ ^2 M3 A* 



SC P 


2 






1.818250 (2) 


1.575289 (2) 


0.701595 (2) 


6.787079 


3 






1.363938 (2) 


1.245279 (2) 




5.218433 


4 






1.235701 (2) 


1.156156 (2) 




4.783715 


5 






1.174520 (2) 


1.114540 (2) 




4.578120 


6 






1.138614 (2) 


1.090432 (2) 




4.458093 


7 






1.114982 (2) 


1.074705 (2) 




4.379373 


SC V 


2 


a = 


45.15762° 


1.606923 


1.504688 


1.493106 (2) 


6.097823 


3 


a = 


45.09672° 


1.280163 


1.228075 


1.223717 (2) 


4.955672 


4 


a = 


45.06976° 


1.182908 


1.147576 


1.145087 (2) 


4.620659 


5 


a = 


45.05455° 


1.135887 


1.109091 


1.107406 (2) 


4.459790 


6 


a = 


45.04479° 


1.108127 


1.086528 


1.085274 (2) 


4.365202 


7 


a = 


45.03799° 


1.089794 


1.071698 


1.070708 (2) 


4.302909 


TABLE V. Parameters of SCs and Wannier indices for 2 + 2e system. The numbers in 
parentheses indicate degree of unstable mode degeneracy. 


Z 


m 


a 


Ml 






A* 


1 


1 


a = 


45° 1.29366 


0.90584 (2) 




3.10533 


2 


1 


a = 


32.2093° 1.36762 


1.33643 




2.70405 


1 


2 


a = 


35.9490° 1.56958 


1.31788 


0.54315 


3.43062 


1 


1836 


a = 


30.0049° 50.32979 


37.46232 




87.79211 
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